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FOREWORD 

This  report  documents  the  development  of  of  a  number  of  techniques  for  the  computer  aided 
analysis  of  quasi-optical  power  combining  systems.  The  major  results  of  the  project  were  the 
development  of  models  for  IMPATT  diodes  suitable  for  harmonic  balance  analysis  and  transient 
analysis;  the  development  of  quasioptical  Green’s  fimctions  for  cavity  resonators  and  grid  systems; 
the  development  of  a  Method  of  Moments  code  which  can  use  the  mixed  scalar  potential  and  vector 
Green’s  fimction  developed  for  quasi-optical  systems;  a  transient  analysis  engine  for  quasi-optical 
systems  and  developments  of  an  harmonic  balance  technique  suitable  for  quasi-optical  systems. 
The  work  makes  a  broader  contribution  to  the  computer  aided  engineering  of  microwave  and 
millimeter  wave  systems  with  complex  device-field  interaictions. 
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1  statement  of  Problem  Studied 

In  quasi-optic2J  power  combining  systems  a  spatial  distribution  of  amplifiers  or  oscillators  and 
radiating  elements  are  used  to  produce  large  amounts  of  power.  The  spatial  distribution  of  the 
amplifiers  and  their  radiating  elements  is  such  that  the  amplifiers  are  distributed  by  a  substantial 
fraction  of  a  wavelength.  The  power  from  the  radiating  elements  is  combined  in  free-space  over  a 
distance  of  many  wavelengths  to  channel  power  into  a  single  paraxial  mode.  The  systems  used  as 
models  in  this  work  are  shown  in  Figures  1  and  2.  Quasioptical  amplifier  and  oscillator  operation 
involves  complex  device-field  interaction.  Computer  aided  simulation  and  subsequent  computer 
aided  design  (CAD)  requires  accurate  modeling  of  the  field  structure  and  a  strategy  to  integrate 
the  field  characterization  into  a  nonlinear  circuit  simulator.  The  nonlinear  simulator  itself  must 
be  capable  of  handling  perhaps  hundreds  of  nonlinear  devices. 

Current  commercial  microwave  computer  aided  engineering  (CAE)  packages  are  oriented  to¬ 
wards  the  design  and  analysis  of  planar  microwave  circuits  including  hybrid  circuits  and  MMIC’s. 
The  computer  aided  analysis  of  these  circuits  is  characterized  by  there  being  relatively  few  active 
devices  cind  by  the  planar  circuit  configuration.  In  nearly  all  situations  the  fields  are  considered 
as  having  a  secondary  effect  to  that  which  can  be  modeled  by  circuit  elements  alone.  Full  field 
simulations  are  also  used  but  the  problem  is  simplified  by  the  circuitry  of  interest  being  confined 
to  one  or  two  planes.  Conventional  microstrip  Greens  functions  can  be  used  in  an  efficient  2^ 
method  of  moments-based  field  solver.  In  contrast,  in  a  quasi-optical  system  the  interaction  of 
fields  with  active  devices  and  with  planar  circuits  is  inherent  to  system  operation.  Furthermore 
quasi-optical  systems  rely  on  combining  the  powers  of  numerous  closely  coupled  active  devi<^. 
Full  field  simulators  have  been  used  to  analyze  individual  cells  of  a  quasioptical  system.  Analyzing 
a  complete  quasioptical  system  using  a  3D  simulator  is  prohibitive  in  terms  of  memory  and  time 
as  a  restdt  of  the  volume  discretization  required  over  many  tens  of  wavelengths.  The  CAE  of  qua¬ 
sioptical  systems  requires  modeling  and  an2dysis  techniques  which  are  not  supported  by  existing 
commercial  microwave  CAE  programs.  The  purpose  of  the  work  proposed  here  is  to  develop  such 
a  modeling  and  analysis  approach  to  support  the  design  and  development  of  quasi-optical  systems. 
The  fundamental  issues  here  are 

1.  handling  device-field  interactions  in  a  non-planar  environment, 

2.  handling  a  very  large  number  of  active  devices  in  steady-state  harmonic  balance  analysis. 

3.  handling  distributed  high  Q  passive  components  in  transient  analysis. 

The  low  efficiencies  achieved  with  current  quasioptical  system  design  methodologies  is,  pre¬ 
sumably,  because  of  poor  impedance  matching  of  the  device  to  the  circuit  in  which  it  is  impeded 
and  perhaps  interaction  with  surrounding  active  devices.  Better  CAE  tools  are  required.  The  aim 
of  this  project  was  to  develop  these  required  CAE  technologies  for  quasi-optical  systems.  Closely 
allied  with  this  was  the  development  of  verification  systems.  An  open  cavity  oscillator,  shown  in 
Figure  1,  and  a  two  dimensional  dielectric  quasioptical  system,  shown  in  Figure  3,  were  developed. 

2  Summary  of  Most  Important  Results 

Contributions  were  made  to  the  computer  aided  engineering  of  quasi-optical  systems  that  were 
realized  in  this  project  aure 
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Figure  1:  Cross-section  of  the  plano-concave  open  resonator. 
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Figure  2:  Grid  amplifier  with  X  and  Y  polarizers. 
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Figure  3:  Grid  amplifier  with  X  aind  Y  polaurizers. 


1.  Development  of  quasi-optical  specific  Green’s  functions.  See  section  2.1. 

2.  Development  of  a  method  of  moments  field  analysis  code  which  uses  the  quasi-optical  Green’s 
functions.  See  section  2.2. 

3.  Development  of  an  X-band  open  cavity  quasi-optical  power  combining  oscillator.  See  section 
2.3. 

4.  Development  of  a  steady-state  (harmonic  balance)  analysis  engine  for  quasi-optical  systems. 
See  section  2.4. 

5.  Development  of  a  transient  analysis  engine  for  high  Q  quasi-optical  systems.  See  section  2.5. 

6.  Development  of  an  enhanced  model  for  and  measurement  based  modeling  technique  for  an 
IMPATT  diode.  See  section  2.6. 

2.1  Quasi-optical  Green’s  Function  Development 

The  approach  that  we  used  considered  the  total  field  radiating  from  a  device  as  a  combination 
of  paraxial  fields  and  nonparaxial  fields.  The  paraxial  fields  are  quasi-optical  fields  described  by 
Hermite- Gaussian  or  Laguerre-Gaussian  modes.  The  nonparaxial  fields  are  described  by  a  free 
space  or  perhaps  half  space  Green’s  function.  In  more  comphcated  situations  with  substrates  the 
free  space  or  half  space  Green’s  function  could  be  replaced  by  the  appropriate  layered  Green’s 
function.  The  trick  is  that  the  paraxial  “free-space  fields”  are  counted  twice.  With  the  double 
accounting  removed  the  complete  field  description  is  obtained.  The  quasioptical  component  of 
the  fields  describes  global  coupling  through  the  cavity  and  the  nonparaxial  fields  describe  near 
neighbor  coupling  via  direct  radiation. 

Details  of  this  work  are  given  in  the  following  publications: 

•  P.L.  Heron,  J.W.  Mink,  G.P.  Monahan,  F.W.  Schwering  and  M.B.  Steer,  “Impedance  ma¬ 
trix  of  an  anteima  array  in  a  quasi-optical  resonator,”  IEEE  Trans.  Microwave  Theory 
Techniques,  Vol.  41,  October  1993,  pp.  1816-1826. 

The  power  from  numerous  nulhmeter  wave  solid-state  sources  can  be  eflSciently  combined 
using  quasi-optical  techniques.  One  such  technique  is  to  place  an  array  of  active  radiating 
sources  within  a  quasi-optical  resonator.  The  driving  point  impedance  of  each  antenna 
is  strongly  effected  by  the  presence  of  all  other  active  antennas  as  well  as  by  the  mode 
structure  and  Q  of  the  resonator.  In  this  paper  the  impedance  matrix  for  an  array  of 
antennas  radiating  into  a  plano-concave  open  resonator  is  determined  through  use  of  the 
Lorentz  integral.  The  resulting  expressions  include  the  effect  of  diffraction  loss  and  are  valid 
for  arbitrary  reflector  spacing,  source  frequency,  array  location  and  geometry.  The  result  can 
be  used  to  impedance  match  each  active  source  to  its  antenna  and  thus  facilitate  design  of  an 
efficient  power  combining  system.  Simulations  using  the  impedance  matrix  in  conjunction 
with  antenna  impedance  model  are  compared  with  two-port  measurements. 

•  P.L.  Heron,  F.W.  Schwering,  G.P.  Monahan,  J.W.  Mink,  and  M.B.  Steer,  “A  dyadic  Green’s 
function  for  the  plano-concave  quasi-optical  resonator,”  IEEE  Microwave  and  Guided  Wave 
Letters,  Vol.  3,  August  1993,  pp.  256-258. 
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Quasioptical  power  combining  using  an  open  resonator  is  an  efficient  and  robust  means  to 
combine  the  power  of  numerous  oscillators.  Here  an  approximate  dyadic  Green’s  function  is 
derived  for  a  quasi-optical  resonator.  The  resonator  consists  of  an  infinite  perfectly  conduct¬ 
ing  planar  reflector  facing  a  partially  reflecting  concave  sphericad  surfaice  of  finite  aperture, 
and  includes  the  effects  of  direct  radiation  loss,  mode  diffraction  losses,  and  conductor  losses. 
Near  field  effects  are  included  so  that  direct  radiation  coupling  between  antennas  can  be  ac¬ 
curately  determined. 

•  Pat  Heron,  Design  of  Millimeter  Wave  Quasi-Optic  Power  Combiners,  Ph.D.  Dissertation, 
North  Cairolina  State  University,  1993. 

Sponsored  research  wais  directed  toward  developing  millimeter  wave  power  sources  utilizing 
quasi-optical  techniques.  A  system  consisting  of  an  airray  of  oscillators  that  radiated  into  a 
quasi-opticail  resonator  wm  analyzed.  Each  oscillator  was  comprised  of  a  solid  state  device 
and  a  radiating  structure.  A  dyaidic  Green’s  function  was  developed  for  a  Fabry-Perot 
resonator  which  consisted  of  a  metallic  planar  reflector  and  a  shallow  spherical  metallic 
reflector.  The  Green’s  function  was  applied  to  determine  the  driving  point  impedance  matrix 
for  an  array  of  electrically  small  antennas  within  the  resonator.  An  experimental  X-band 
resonator  was  designed  and  fabricated,  then  one  and  two-port  measurements  were  used  to 
validate  the  theoretical  calculations.  A  technique  wais  determined  for  simulation  of  antennas 
that  «ire  not  electrically  smadl  which  radiate  into  the  cavity.  These  techniques  are  shown  to 
be  applicable  to  coupling  structures  for  quasi-optical  systems  in  general. 

Practical  considerations  regairding  the  simulation  of  nonlinear  solid  state  driving  elements 
were  addressed.  A  a  technique  for  efficient  Jacobian  calculation  using  the  multidimensionaJ 
fast  Fourier  transform  as  well  as  a  technique  for  simulator  time-domain  oversampling  were 
developed  so  that  multiple  oscillator  systems  can  be  efficiently  simulated.  These  simulation 
techniques  were  implemented  and  tested  in  a  harmonic  balance  circuit  simulator. 

•  P.L.  Heron,  G.P.  Monahan,  F.W.  Schwering,  J.W.  Mink,  and  M.B.  Steer  “Multiport  circuit 
model  of  «in  antenna  array  in  an  open  quasi-optical  resonator,”  Proc.  1993  URSI  Conference, 
Jime  1993,  p.  84. 

In  a  quasi-optical  power  combiner  the  output  of  many  solid-state  sources  is  efficiently  aind 
robustly  combined.  This  is  cichieved  «is  a  resffit  of  the  one-to-many  coupling  as  the  power 
radiated  by  a  single  element  is  reflected  by  the  curved  reflector  on  to  many  oscillators. 
The  driving  point  impedance  of  each  eintenna  is  strongly  effected  by  the  presence  of  all  other 
cictive  antennas  as  well  as  by  the  mode  structure  and  Q  of  the  resonator.  In  this  presentation 
the  derivation  of  a  field  theoretic  model  describing  the  coupling  is  described.  An  impedance 
matrix  is  determined  through  use  of  the  Lorentz  reciprocity  theorem  and  cadculation  of  the 
aperture  and  conductor  losses.  A  half  space  Green’s  function  is  used  to  incorporate  the  near¬ 
field  effect  of  non  quasi-optical  fields.  The  result  obtauned  can  be  used  to  impedaince  match 
each  active  source  to  its  antenna  and  thus  facilitate  design  of  aui  efficient  power  combining 
system.  As  am  exaunple  the  cadculated  and  meaisured  transimpedamce  of  two  amtennas  are 
considered. 

•  P.L.  Heron,  G.P.  Monaham,  J.E.  Byrd,  M.B.  Steer,  F.W.  Schwering  and  J.W.  Mink,  “Circuit 
level  modeling  of  quaisiopticaJ  power  combining  open  cavities,”  1993  IEEE  MTT-S  Interna¬ 
tional  Microwave  Symposium  Digest,  June  1993,  pp.  433-436. 
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A  multiport  circuit  level  model  is  developed  for  an  open-cavity,  quasioptical  power  combiner. 
The  model  is  developed  using  Hermite- Gaussian  beam  mode  theory,  the  Lorentz  reciprocity 
theorem  and  determination  of  diffraction  losses.  An  impulse  response/ convolution-based 
technique  being  developed  for  incorporating  this  high  Q  model  in  a  transient  simulator  is 
discussed. 

2.2  Method  of  Moments  Analysis 

All  other  schemes  for  hamdling  the  fields  in  quasi-optical  systems  use  a  imit  ceU  approach  where 
ideal  magnetic  and  electric  current  walls  are  used.  The  unit  cell  based  modeling  approach  is  akin 
to  the  circuit  element  modeling  in  conventional  microwave  simulators.  While  this  leads  to  fast 
circuit  simulation  compromises  are  made.  This  Green’s  function  uses  a  mixed  paraxial/nonparaxial 
Green’s  function  approach  which  is  combined  with  an  extremely  efficient  Conventional  Method 
of  Moments  code  using  uses  a  scalar  Green’s  function  whereas  the  quasioptical  Green’s  function 
is  a  mixed  scalar  potential  and  vector  Green’s  function.  It  would  be  prohibitively  complex  to 
develop  a  scalar  potential  Green’s  function  for  a  quasi-optical  system  including  effects  such  as 
diffraction  losses  and  conductive  losses.  However  the  technique  that  we  developed  supports  the 
heiiristic  development  of  Green’s  functions  and  not  just  for  quasi-optical  systems.  This  could  lead 
to  important  developments  in  the  top-down  design  of  microwave  and  millimeter- wave  systems  in 
that  an  approximate,  heuristically-based.  Green’s  function  could  be  used  in  early  system  design. 

Details  of  this  work  are  given  in  the  following  publications: 

•  T.W.  Nuteson,  G.P.  Monahan,  M.B.  Steer,  K.  Naishadham,  J.W.Mink,  K.  Kojucharoff  and 
J.  Harvey,  “Full-wave  analysis  of  quasi-optical  structures,”  IEEE  Trans.  Microwave  Theory 
Techniques,  1996,  In  Press. 

A  full-wave  moment  method  implementation,  using  a  combination  of  spatial  and  spectral 
domains,  is  developed  for  the  analysis  of  quasi-optical  systems.  An  electric  field  dyadic 
Green’s  function,  including  resonant  and  non-resonant  terms  corresponding  to  coupling  from 
modal  and  non-modal  fields,  is  employed  in  a  Galerkin  routine.  The  dyadic  Green’s  function 
is  derived  by  separately  considering  paraxial  and  non-paraxi£d  fields  and  is  much  easier  to 
develop  than  a  mixed,  scalar  and  vector,  potential  Green’s  function.  The  driving  point 
impedance  of  several  antenna  elements  in  a  quasi-optical  open  cavity  resonator  and  a  3  x  3 
grid  in  free  space  are  computed  and  compared  with  measurements. 

•  T.W.  Nuteson,  G.P.  Monahan,  M.B.  Steer,  K.  Naishadham,  J.W.  Mink,  and  F.K.  Schwering, 
“Use  of  the  moment  method  and  dyadic  Green’s  functions  in  the  analysis  of  quasi-optical 
structures,”  1995  IEEE  MTT-S  International  Microwave  Symposium  Digest,  May  1995.  pp. 
913-916. 

A  moment  method  using  a  dyadic  Green’s  function  is  developed  for  the  analysis  of  quasi¬ 
optical  systems.  The  dyadic  Green’s  function  used  has  separate  terms  for  the  paraxial  and 
non-paraxial  fields  and  is  much  easier  to  develop  than  a  mixed  potential  Green  s  function. 
The  method  is  applied  to  the  analysis  of  antenna  elements  in  a  quasi-optical  resonator. 


2.3  X-band  Open  Cavity  Quasi-Optical  Power  Combining  Oscillator 

An  X-band  open  cavity  quasi-optical  oscillating  power  combining  system,  as  shown  in  Figure  1, 
vfas  constructed  as  one  of  the  verification  systems.  The  oscillator  used  IMPATT  diodes  but  the 
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results  were  disappointing. 

Detfiils  of  the  work  axe  described  in 

•  Gregory  Monahan,  Characterization  and  Design  of  Millimeter  Wave  Quasi-Optical  Power 
Combiners,  Ph.D.  Dissertation,  North  Carolina  State  University,  1995. 

Research  was  directed  towards  developing  an  experimentally  based  design  methodology  for 
quasi-optical  power  combining  using  IMPATT  diode  oscillators.  The  design  and  character¬ 
ization  of  an  open  cavity  resonator  for  use  in  a  quasi-optical  power  combiner  is  presented. 
The  lowest  order  Hermite  Gaussian  and  Laguerre  Gaussian  modes  of  Fabry-Perot  (open 
cavity)  resonators  are  shown  to  be  degenerate.  A  coaxial  test  fixture  was  designed,  fabri¬ 
cated  and  used  to  experimentally  characterize  post  mounted  IMPATT  diodes  in  the  small 
signal  regime.  This  i^ormation  was  then  used  to  design  a  radiating  oscillator  operating  at 
8.65  GHz  consisting  of  an  IMPATT  diode  connected  directly  to  a  rectangular  patch  antenna. 
The  operation  of  the  oscillator  in  free  space  and  in  the  presence  of  cavity  modes  was  charac¬ 
terized.  The  effect  of  the  cavity  modes  upon  the  operation  of  the  oscillator  was  explored.  An 
array  of  three  IMPATT  diode  oscillators  operating  in  the  open  cavity  resonators  was  devel¬ 
oped  using  the  design  methodology  presented  herein.  Design,  construction,  instrumentation, 
measurements,  and  data  presentation  of  the  experimental  apparatus  is  described. 

•  G.P.  Monahan,  P.L.  Heron,  M.B.  Steer,  J.W.  Mink  and  F.W.  Schwering  “Mode  degeneracy 
in  quasi-optical  resonators,”  Microwave  and  Optical  Technology  Letters,  Vol.  8,  No.  5,  April 
5  1995,  pp.  230-232. 

Fabry-Perot  resonators  cind  beam  waveguides  are  used  in  quasi-optical  systems  to  refocus 
electromagnetic  fields.  These  fields  are  described  by  Laguerre  Gaussian  or  Hermite  Gaussian 
quasi-TEM  mode  families  dependent  on  whether  the  aperture  is  rectangular  or  circular.  It  is 
shown  that  certain  lower  order  Laguerre  Gaussian  and  Hermite  Gaussian  modes  are  identical. 
This  has  implications  for  the  design  of  quasi-optical  systems  for  mode  selection. 

2.4  Steady-State  (Harmonic  Balance)  Analysis  of  Quasi- Optical  Power 
Combiners 

Modifications  were  made  to  the  conventional  harmonic  balance  algorithm  to  support  active  quasi- 
optical  systems.  The  2uialysis  engine  incorporates  novel  memory  and  anailysis  efficiencies  so  that 
the  ^lnaly3is  of  lOO’s  of  active  devices  in  a  quasioptical  system  is  possible.  Procedures  for  the 
optimization  of  power  added  efficiency  and  output  power,  subject  to  stable  optimization,  will  be 
developed.  Currently  all  of  the  elements  must  be  assumed  to  oscillate  at  the  same  frequency. 
The  CAD  tool  enables  the  following  to  be  investigated:  the  phase  distribution  across  the  antenna 
array;  power  tind  oscillation  frequency. 

Background  for  the  work  is  given  in  the  following  publications: 

•  M.B.  Steer  and  S.G.  Skaggs,  “CAD  of  GaAs  microwave  circuits:  historical  perspective  and 
future  trends,”  invited  paper,  Korea-United  States  Design  and  Manufacturing  Workshop, 
November  17-19,  1993,  Taejon,  Korea,  pp.  109-118. 

The  nonlinear  analysis  of  microwave  circuits  has  reached  a  high  level  of  maturity  over  the  last 
decade.  By  assuming  that  only  a  fimte  number  sinusoids  are  present  in  a  nonlinear  circuit, 
the  computational  burden  of  computing  the  transient  response  of  the  circuit  is  avoided  and 
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only  the  steady  state  response,  given  by  the  amplitudes  and  phases  of  the  sinusoids,  is 
required.  The  Harmonic  Balance  technique  is  ideally  suited  to  solving  for  the  phasors  of 
the  steady  state  response.  However  this  technique  is  not  as  versatile  as  transient  analysis 
techniques.  An  historical  perspective  is  presented  and  future  trends  and  requirements  for 
the  CAD  of  GaAs  microwave  circuits  are  outlined. 

•  M.B.  Steer  and  S.G.  Skaggs,  “CAD  of  GaAs  microwave  circuits:  historical  perspective  and 
future  trends,”  invited  paper,  Korea-United  States  Design  and  Manufacturing  Workshop, 
November  17-19,  1993,  Taejon,  Korea,  pp.  109-118. 

Details  of  this  work  that  have  not  been  published  elsewhere  are  given  in  Appendix  A. 


2.5  Transient  Analysis  of  a  Quasi-Optical  Power  Combiner 

Current  design  and  operation  concerns  of  quasi-optical  combining  systems  largely  involve  locking 
behavior.  Much  of  the  concept  behind  quasi-optical  power  combining  relies  on  global  locking 
of  the  various  oscillating  elements.  However  maximizing  the  power  extracted  from  the  system 
requires  that  the  elements  be  placed  relatively  closely  to  each  other.  This  increases  the  nearest 
element  interactions  so  that  local  locking  rather  than  global  locking  becomes  dominant.  However 
some  local  locking  is  required  to  ensure  that  the  sources  initially  lock.  To  model  these  effects 
transient  simulation  is  required.  We  have  used  an  existing  convolution-based  transient  simulator 
to  successfully  model  a  simple  model  of  a  single  element  quasi-optical  cavity  oscillator.  The 
simulation  strategy  uses  the  impulse  response  of  the  linear  circuitry  to  interface  the  frequency- 
dependent  cavity  model  to  the  nonlinear  oscillating  element  model.  Interfacing  is  achieved  through 
convolution. 

Transient  analysis  is  essential  in  determining  some  of  the  most  critical  parameters  of  a  quasiop- 
tical  systems.  Locking  behavior  and  O^occurrence  of  chaotic  behavior  cannot  be  determined  any 
other  way.  However  transient  analysis  of  steady  state  microwave  systems  is  notoriously  slow.  Far 
too  slow  for  design  optimization,  hi  conventional  microwave  circuit  design  the  lengthy  simulation 
time  is  reduced  through  application  of  harmonic  balance  analysis.  Simulation  time  reduction  by 
a  factor  of  a  thousand  (and  very  often  much  more)  is  obtained.  It  then  becomes  feasible  to  apply 
circuit  optimization  techniques  to  maximize  power  jidded  efficiency,  output  power  and  other  design 
parameters. 

We  have  adapted  a  technique  developed  in  another  project  (M.S.  Basel,  M.B.  Steer  and  P.D. 
Franzon  “Simulation  of  High  Speed  Interconnects  Using  a  Hierarchical  Packaging  Simulator,” 
IEEE  Trans,  on  Components  Hybrids  and  Manufacturing  Technology/Advanced  Packaging,  Febru¬ 
ary  1995,  pp.  74-82.)  for  the  transient  analysis  of  distributed  circuits  and  in  particular  for  high 
Q  quasi-optical  systems.  The  following  describes  what  was  done:  The  transient  analysis  of  a 
distributed  circuit  without  a  canonical  model  requires  the  use  of  Fourier  transform  techmques  to 
translate  the  frequency  domain  charaicterization  of  the  linear  circuit  to  an  impulse  response  that 
can  be  used  in  a  transient  convolution-based  analysis.  This  transformation  would  normally  lead 
to  significant  aliasing  for  high  Q  circuits.  In  the  work  by  Basel  et  al.  we  presented  a  method 
which  preserves  the  characterization  in  generating  a  modified  impulse  response.  The  modification 
of  the  impulse  response  is  removed  during  iteration.  Prototype  code  for  performing  the  analysis 
of  a  quasioptical  system  has  been  completed. 

In  transient  analysis  either  a  canonical  model  of  the  distributed  quasioptical  system  must  be 
developed  or  else  an  impulse  response  derived  by  Fourier  transforming  a  very  large  number  of 
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model  evaluations  at  discrete  frequencies.  Modeling  compromises  are  made  in  both  cases.  The 
canonical  modeling  approach  models  the  complicated  distributed  circuit  (the  quasioptical  system) 
by  a  lumped  element  model  In  the  impulse  response  approach  which  we  use  the  aliaising  effects  can 
never  be  entirely  eliminated  and  many  more  model  evaluations  will  invariably  be  required  than 
that  required  in  an  harmonic  balance  analysis. 

Details  of  this  work  are  mostly  given  in  the  paper  by  Basel  et  al.  Prototype  code  has  been 
written  but  the  code  is  not  compatible  with  the  Spice  simulation  approeich.  This  will  be  addressed 
in  the  subsequent  commercialization  activities. 

Backgroimd  emd  description  of  the  work  conducted  early  in  this  project  are  described  in  the 
following  publication: 

•  M.B.  Steer,  “Simulation  of  microwave  and  millimeter- wave  oscillators,  present  capability  and 
future  directions,”  Proc.  1992  Workshop  on  Integrated  Nonline£ir  Microwave  and  Millimeter- 
wave  Circuits,  October  1992  (invited  keynote  paper). 

The  current  state-of-the  art  of  oscillator  simulation  techniques  is  presented.  C<indidate 
appro2w:hes  for  the  next  generation  of  oscillator  simulation  techniques  are  reviewed.  The 
method  is  presented  which  uses  an  efficient  and  robust  convolution-based  procedure  to  in¬ 
tegrate  frequency- domain  modeling  of  a  distributed  linear  network  in  trmsient  simulation. 
The  impulse  response  of  the  entire  linear  distributed  network  is  obtained  and  the  algorithm 
presented  herein  ensures  that  aliasing  effects  are  minimized  by  introducing  a  procedure  that 
ensures  that  the  interconnect  network  response  is  both  time-limited  eind  bcind-limited.  In 
particular,  artificial  filtering  to  bandlimit  the  response  is  not  required. 

Details  of  the  work  that  have  not  been  published  elsewhere  are  given  in  Appendix  B. 

2.6  IMPATT  Device  Model 

a  The  open  cavity  power  combining  oscillator  used  IMPATT  diodes  as  the  active  device  and  it 
was  anticipated  that  an  existing  IMPATT  diode  model  could  be  used  in  computer  simulations. 
However  the  existing  models  do  not  conserve  charge  when  they  are  operated  in  the  large  signal 
regime.  As  such  the  development  of  Iwge  signed  IMPATT  diode  models  suitable  for  transient 
analysis  and  heirmonic  balemce  analysis  were  developed.  Extensive  eictive  device  modeling  was  not 
anticipated  in  this  project. 

Details  of  this  work  that  have  not  been  published  elsewhere,  aue  given  in  the  Appendix  C 

Experimental  IMPATT  Characterization 

Detadls  of  the  experimental  IMPATT  characterization  work  are  given  in  the  following  publica¬ 
tions: 

•  M.B.  Steer,  “Diode  characterization  in  a  coaxiad  mount,”  International  Journal  of  Microwave 
and  Millimeter  Wave  Computer  Aided  Engineering,  Vol.  3,  April  1993,  pp.  114-117. 

A  previously  presented  technique  for  experimentadly  characterizing  packaged  microwave 
diodes  mounted  in  a  coaxial  test  fixture  is  investigated.  It  is  demonstrated  that  the  fields 
surrounding  the  diode  and  the  cadibration  packages  are  awiequately  developed  to  provide  am 
accurate  diode  equivadent  circuit  model  for  use  computer  added  design. 

•  G.P.  Monahain,  A.S.  Morris  amd  M.B.  Steer,  “A  coaodad  test  fixture  for  characterizing  low 
impedance  microwave  two-terminad  devices,”  Microwave  and  Optical  Technology  Letters,  Vol. 
6,  Mairch  1993,  pp.  197-200. 
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A  coaxial  test  fixture  is  designed  and  used  to  measure  the  input  impedance  of  a  packaged 
IMP  ATT  diode  in  the  frequency  ramge  from  0  to  18  GHz.  In  this  paper  such  a  mount 
is  presented  which  incorporates  a  coaxiad  transformation  to  transform  the  relatively  high 
measurement  characteristic  impedaince  to  low  impedance  of  a  coaxial  line  with  atn  oversized 
center  conductor.  Following  the  coaxial  transformer  the  fields  aire  converted  to  a  raulial  struc¬ 
ture  corresponding  to  the  fields  of  a  mounted  diode.  This  diode  characterization  procedure 
is  being  used  to  develop  IMP  ATT  diode  models  for  use  in  the  simulation  of  a  quasioptical 
power  combining  oscillator  wherein  the  diode  appears  under  a  sub-resonant  patch  antenna. 
The  work  is  an  extension  of  our  earlier  work  but  with  greater  attention  paid  to  the  coaxial 
transformer. 

2.7  Two-Dimensional  Quasi-Optical  Power  Combiner 

•  H.  Hwang,  T.W.  Nuteson,  M.B.  Steer,  J.W.  Mink,  J.  Harvey  and  A.  Paolella,  “Quasi-optical 
power  combining  in  a  dielectric  substrate.  International  Symposiiim  on  Signals,  Systems  and 
Electronics,  October  25-27,  1995,  pp.  89-92. 

This  paper  presents  a  quasi-optical  MESFET-based  dielectric  slab  waveguide  with  an  am¬ 
plifier  array.  Up  to  19  dB  amplifier  gain  and  9.5  dB  system  gain  is  obtained  at  7.38  GHz. 
Measurements  of  amplifier  and  system  gain,  output  power  vs.  input  power  and  trans¬ 
verse  power  distribution  are  presented. 

•  H.  Hwang,  G.P.  Monahan,  M.B.  Steer,  J.W.  Mink,  and  F.K.  Schwering,  “A  dielectric  slab 
waveguide  with  four  planar  power  amplifiers,”  1995  IEEE  MTT-S  International  Microwave 
Symposium  Digest,  May  1995.  pp.  921-924 

A  hybrid  dielectric  slab  beam  waveguide  with  four  MESFET  amplifiers  employing  quasi- 
optical  power  combining  is  reported  for  the  first  time.  Up  to  10  dB  power  gain  is  obtained 
at  7.4GHz.  Measurements  for  power  gain,  amplifier  gain,  insertion  loss  and  transverse  power 
distribution  are  presented.  The  fabrication  technique  employed  is  suitable  for  planar  MMIC 
circuits. 

•  F.  Poegel,  S.  Irrgang,  S.  Zeisberg,  A.  Schuenemann,  G.P.  Monahan,  H.  Hwang,  M.B.  Steer, 
J.W.  Mink,  and  F.  K.  Schwering,  “Demonstration  of  an  oscillating  quasi-optical  slab  power 
combiner,”  1995  IEEE  MTT-S  International  Microwave  Symposium  Digest,  May  1995.  pp. 
917-920. 

Power  combining  in  a  hybrid  dielectric  slab-beam  waveguide  resonator  using  a  MESFET 
osciUator  array  is  reported  for  the  first  time.  Four  MESFET  osciUators  lock  via  quasi-optical 
modes  to  produce  a  signal  at  7.4  GHz  with  a  3  dB  linewidth  of  less  than  3  kHz. 

•  A.  Schuneman,  S.  Zeisberg,  P.  L.  Heron,  G.  P.  Monahan,  M.  B.  Steer,  J.  W.  Mink,  F.  W.  Schw¬ 
ering,  “A  prototype  quasi-optical  slab  resonator  for  low  cost  millimeter-wave  power  com¬ 
bining,”  Workshop  on  Millimeter- Wave  Power  Generation  and  Beam  Control,  Huntsville, 
Alabama,  September  14-  16,  1993,  pp.  235-243. 

A  quasi-optical  slab  resonator  is  a  planar  version  of  an  open  quasi-optical  resonator  and  has 
many  potential  manufacturabiUty  and  cost  advantages.  It  uses  a  waveguiding  mechanism 
that  combines  beam- waveguide  and  dielectric  waveguide  propagation  modes  to  confine  fields 
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to  a  plane.  Here  the  slab  resonator  is  experimentally  characterized  and  design  data  presented 
for  a  prototype  quasi-optical  resonator. 

•  J.W.  Mink,  F.W.  Schwering,  P.L.  Heron,  G.P.  Monahan,  A.  Schuneman,  S.  Zeisberg,  and 
M.B.  Steer,  “A  hybrid  dielectric  slab-beam  waveguide,  theory  and  experiment,”  19  th  Army 
Science  Conference^  June  1994. 

•  A.  Schuneman,  S.  Zeisberg,  G.P.  Monadian,  F.W.  Schwering,  J.W.  Mink,  and  M.B.  Stwr, 
“Experimental  investigation  of  a  quasioptical  slab  resonator,”  Proc.  1993  URSI  Conference, 
June  1993,  p.  27. 

Quasi-optical  power  combiners  using  open-resonators  achieve  efficient  and  robust  combining 
at  millimeter- wave  frequencies  and  above.  However  systems  cannot  be  photolithographically 
defined  thus  limiting  mass  production  and  contributing  to  high  cost.  A  similar  structure  more 
amenable  to  photolithographic  reproduction  is  the  hybrid  dielectric  slab-beam  waveguide 
(HDSBW).  This  structure  combines  the  waveguiding  principles  of  dielectric  surface  waves  of 
a  slab  guide  surface  and  the  confined  beam  corresponding  to  Gauss- Hermite  beam  modes. 
In  this  paper  a  HDSBW  slab  resonator  was  experimentally  investigated. 

•  S.  Zeisberg,  A.  Schunemann,  G.P.  Monahan,  P.L.  Heron,  M.B.  Steer,  J.W.  Mink  and  F.W. 
Schwering,  “Experimental  investigation  of  a  quasi-optical  slab  resonator,”  IEEE  Microwave 
and  Guided  Wave  Letters,  Vol.  3,  August  1993,  pp.  253-255. 

A  quasi-optical  slab  resonator  for  TE  modes  was  experimentally  charaicterized  to  demon¬ 
strate  a  plwar  technology  for  quasi-optical  devices.  Predicted  and  measured  frequencies  of 
resonance  of  the  TE  slab  modes  and  electric  field  profiles  are  in  close  agreement. 

•  H.  Hwang,  T.W.  Nuteson,  M.B.  Steer,  J.W.  Mink,  J.  Harvey  and  A.  PaoleUa,  “Slab-based 
quasi-optical  power  combining  system,”  Infraired  and  Millimeter- Wave  Conference,  Orlando 
Florida,  December  1995. 

A  slab-b^ed  quasi-optical  power  combining  system  with  convex  and  concave  lenses  is  inves¬ 
tigated.  Experimental  results  imply  that  a  concave-lens  system  has  less  scattering  loss  and 
higher  system  gains  than  a  concave-lens  system.  An  amplifier  gain  of  15  dB  and  a  system 
gain  of  8  dB  were  achieved. 

•  H.  Hwang,  T.W.  Nuteson,  M.B.  Steer,  J.W.  Mink,  J.  Harvey  and  A.  P<iolella,  “Quasi-optical 
power  combining  techniques  for  dielectric  substrates,”  International  Semiconductor  Device 
Research  Symposium,  December  6-8,  1995. 

The  dielectric  slab  system  described  here  has  the  advantage  of  being  two-dimensional  eind 
is  thus  more  amenable  to  photolithographic  reproduction  than  the  conventional  open  quasi- 
optical  power  combining  structures.  The  previous  investigations  of  the  quasi-optical  di¬ 
electric  slab  cavity  and  wavegtiide  demonstrated  the  suitability  of  this  structure  for  the 
integration  of  quasi-optical  power  combining  with  MMIC  technology.  A  complete  DSBW 
quasi-optical  system,  as  shown  in  Fig.  1,  could  consist  of  the  following:  a  source,  active  or 
injection,  an  zimplifier  array,  triplers,  cind  a  leaky  wave  antenna  which  would  be  used  for 
steering  the  energy  out  of  the  system.  Between  each  of  the  stages  lenses  are  used  to  focus 
the  guided  waves  for  optimad  field  concentration  on  the  elements  in  the  system.  In  this 
work  we  present  the  amplifier  array  stage  using  both  convex  and  concave  lenses  as  shown 
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in  Fig.  2.  The  DSBW  amplifier  system  incorporates  four  MESFET  amplifiers  and  two  thin 
convex/concave  lenses.  The  waveguide  system  was  adjusted  with  the  transistors  turned  off 
so  that  the  guided  waves  are  focused  near  the  aperture  of  the  receiving  horn.  The  dielectric 
slab  is  Rexolite  (e  =  2.57,  tan^  =  0.0006  at  X-band  ),  and  it  is  27.94  cm  wide,  62  cm  long, 
and  1.27  cm  thick.  The  convex  lenses  are  fabricated  from  Macor  (c  =  5.9,  tan^  =  0.0025 
at  100  kHz)  with  a  radius  of  30.48  cm,  and  the  focal  length,  /,  is  28.54  cm.  The  concave 
lenses  are  air  (e  =  1)  with  a  radius  of  30.48  cm,  and  the  focal  length,  /,  is  40.4  cm.  The 
aperture  width  of  both  horn  antennas  is  9  cm,  designed  to  be  wide  enough  to  catch  most  of 
the  amplified  power.  Energy  emitted  from  the  input  radiator  propagates  in  a  quasi-optical 
TE  Gaussian  mode  in  the  dielectric  slab  waveguide, and  is  focused  by  the  first  lens  in  the 
middle  area  of  the  slab.  This  system  is  designed  so  that  the  amplifier  unit  cells  are  within 
the  beam  waist  (the  1/e  field  points). 
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4  Future  Direction 

The  work  described  here  is  being  continued  in  three  projects: 

1.  “Circuit  Modeling  and  Computer  Aided  Design  for  Quasioptical  Systems,”  Scientific  Re¬ 
search  Associated,  Inc.  under  U.S.  Army  Missile  Command  contract  DAAH01-95-C-R111. 

2.  “Computer  Aided  Engineering  Tools  for  Microwave  and  Millimeter- Wave  Quasi-optical  Os¬ 
cillators  and  Amplifiers,”  Compact  Software,  Inc.  under  ARPA  SBIR  Program  contract 
DAAH01-95-C-R203. 

3.  “Research  and  Development  of  a  Quasi-Optical  Dielectric  Slab  Power  Combining  System,” 
U.S.  Army  Research  Office  through  grant  DAAH04-95-1-0536. 

The  first  two  projects  involve  commercialization  of  the  activities  conducted  under  the  project 
reported  on  here.  The  third  project  is  continuation  of  the  two-dimensional  power  combiner  work 
described  here  and  the  development  of  a  complete  two-dimensional  power  combining  system. 

A  Harmonic  Balance  Analysis  of  Quasi-Optical  Power 
Combining  Circuits 

Harmonic  balance  analysis  has  become  the  accepted  method  for  finding  the  steady  state  response 
of  nonlinear  microwave  circuits.  This  analysis  requires  the  solution  of  a  system  of  nonlinear 
equations  by  Newton’s  method,  and  is  often  limited  by  the  size  of  the  Jacobian  matrix.  For 
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a  large  matrix,  both  memory  usage  and  run  time  costs  can  become  significant.  $  This  papers 
commercial  harmonic  balcince  simulators  use  different  sparsification  techniques.  The  essence  of 
the  work  we  propose  here  is  to  merge  the  sparsification  techniques  to  obtain  a  efficient  harmonic 
balance  simulation  for  circuits  with  meiny  active  devices. 

Harmonic  balance  analysis  is  a  robust  and  efficient  technique  for  the  analysis  of  small  circuits, 
but  memory  demands  increase  rapidly  as  problem  size  increases.  Memory  usage  is  a  critical  factor 
for  very  large  problems  as  once  RAM  is  exceeded,  disk  swapping  is  used  with  an  associated  feictor 
of  10  reduction  in  speed.  In  modern  harmonic  balance  programs,  in  which  the  elements  of  the 
Jacobi£in  Me  calculated  analytically,  the  model  evaluation  time  is  comparable  to  that  required  for 
matrix  formulation  (R.  Gilmore  and  M.B.  Steer,  “Nonlinear  Circuit  Analysis  Using  the  Method  of 
Harmonic  Balance  —  a  Review  of  the  Art:  Part  II,  Advanced  Concepts”  International  Journal  on 
Microwave  and  Millimeter  Wave  Computer  Aided  Engineering,  Vol.  1,  April,  1991,  pp.  159-180) 
For  large  problems,  for  example  with  many  tones  and/or  many  nonlinear  elements,  the  matrix 
manipulation  time  becomes  the  major  component  of  the  simulation  time.  Also,  the  size  of  the 
Jacobian  matrix  dominates  memory  usage.  With  the  aim  of  reducing  both  memory  utilization 
and  simulation  time,  conventional  and  blocked  sparse  matrix  techniques  and  Samanskii’s  (chord) 
method  have  been  used.  Each  approach  to  sparsification  dramatically  reduces  memory  usage  and 
simulation  times.  However  the  memory  usage  and  simulation  times  are  still  extremely  large  for 
many-device  circuits  such  as  Quasioptical  <implifiers. 

Technically,  the  modifications  require  combining  the  block  Newton  iteration  method  (C.R. 
Chang,  P.L.  Heron  and  M.B.  Steer  “Hjurmonic  Balance  and  Frequency  Domain  Simulation  of 
Nonlinear  Microwave  Circuits  using  the  Block  Newton  Method,”  IEEE  Trans.  Microwave  Theory 
Techniques,  April  1990,  pp.  431-434.)  efficient  Jacobi2tn  calculation  (P.L.  Heron,  and  M.B.  Steer, 
“Jacobian  calculation  using  the  multidimensional  feist  Fourier  transform  in  the  harmonic  balance 
analysis  of  nonlinear  microwave  circuits,”  IEEE  Trans.  Microwave  Theory  Techniques,  April  1990, 
pp.  429-431)  with  the  conventional  sparse  matrix  technology.  Extensions  to  detect  instabilities 
leading  to  chaotic  phenomena  in  such  a  large  system  remain  to  be  developed. 

A.l  Newton  Iteration  Procedure 

In  applying  harmonic  balaince  analysis,  a  system  of  nonlinear  circuit  equations  is  developed  and  am 
error  fimction  is  formulated  and  minimized.  This  error  function  typically  involves  the  Kirchhoff’s 
current  law  error  at  the  interface  between  the  linear  cind  nonlinear  circuits,  but  may  include  other 
error  determinants  (such  as  Kirchhoff’s  voltage  law)  for  those  elements  that  can  not  be  described  as 
controlled  current  sources.  The  error  function  is  minimized  with  respect  to  a  set  of  state  variables 
to  an  acceptable  level  using  Newton’s  method.  The  state  variables  are  typically  the  phasors  of 
the  voltages  at  the  interface  nodes,  although  the  use  of  alternative  state  variables  dramatically 
improves  convergence  in  some  cases. 

The  harmonic  balance  problem  can  be  formulated  in  an  iterative  form  as  follows 

=’  X  —  aJ“^('x)  •  f(’x).  (1) 

where  f(*x)  the  error  vector  at  the  iterate  and  *x  the  vector  of  state  variables  at  the  iterate 
of  f(*x),  the  iterative  update  procedure  is  given  by 

The  harmonic  b«ilance  error  is  defined  as  the  norm  of  the  f('x)  vector  in  (1)  and  iteration  continues 
until  this  error  is  less  than  some  user-defined  value. 

Various  modifications  of  the  Newton  update  procedure  outlined  in  (1)  Me  possible.  Rarely  is 
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(2) 


the  inverse  of  the  Jticobian  evaluated  directly.  Instead  the  corresponding  matrix  equation 

J('x)(’'’'^x  — '  x)  =  f(’x) 

is  solved  for  ’"•'^x  using  LU  decomposition.  The  '■*'^x  estimate  is  then  adjusted  by  the  scaling  factor 
a.  From  the  solution  of  (2)  J“^(’x)  can  be  evaluated.  Furthermore,  in  Samanskii’s  modification 
of  the  procedure,  J“^(’x)  determined  at  one  iteration  is  reused  at  subsequent  iterations  as  long  as 
the  harmonic  balance  error  decreases.  In  practice  J(x)  is  a  sparse  matrix  and  many  of  the  entries 
are  relatively  small.  Through  ordering  of  the  rows  and  columns  of  the  Jacobian,  the  significant 
entries  form  blocks  along  the  diagonal  with  only  relatively  small  entries  outside  these  blocks.  By 
ignoring  the  out-of-block  entries  the  Jacobian  becomes  a  block-diagonal  Jacobian  which  can  be 
manipulated  at  a  fraction  of  the  memory  and  time  costs  that  a  full  Jacobian  can  be  manipulated. 
The  sparsity  can  be  taken  advantage  of  using  conventional  sparse  matrix  representation  resulting 
in  dramatic  reductions  in  memory  usage  and  increases  in  computation  speed.  The  performan^ 
improvements  obtained  using  different  combinations  of  sparse  matrix,  blocking,  and  Samanskii’s 
method  are  the  subject  of  this  paper.  Greater  specifics  are  discussed  below  in  the  context  of  the 
specific  implementation. 

A.2  Results  to  Date 

Aspects  of  the  proposed  approach  have  been  implemented  in  a  prototsrpe  simulation  program 
called  TRANSIM.  First,  TRANSIM  forms  the  modified  nodal  admittance  matrix  of  the  circuit, 
which  is  used  to  obtain  the  linear  response  to  each  iterate  of  x.  Next,  the  interface  unknowns  are 
transformed  to  the  time  domain  and  applied  to  the  nonlinear  elements  to  determine  the  nonlinear 
response  and  error  vector.  If  the  resulting  error  is  less  than  the  error  of  the  previous  iterate,  the 
Jacobian  is  calculated.  Otherwise,  the  scaling  factor  a  is  reduced  and  the  same  Jacobian  is  reused. 
The  time  domain  derivatives  are  converted  directly  to  frequency  domain  derivatives  by  use  of  a 
multidimensional  fast  Fourier  transform(NFFT)  technique. 

As  seen  in  (1),  the  term  J'H’x)  ’  f(‘*)  I^^stead  of  performing  the  inversion  of 

J,  TRANSIM  solves  the  Unear  equation  J(’x)  •  b  =  /(*x)  by  LU  factorization  of  J(*x).  Since  the 
LU  factorization  of  a  matrix  A  need  be  computed  only  once  to  solve  any  equation  Ax  =  b,  LU 
factorization  can  help  reduce  run  time  when  the  Jacobian  is  reused. 

Perhaps  the  simplest  way  to  reduce  run  time  is  the  method  proposed  by  Samanskii,  also  called 
the  chord  method.  Since  the  Jacobian  often  does  not  change  significantly  between  iterations,  it  is 
not  always  necessary  to  recalculate  the  Jacobian  at  each  step.  Instead,  the  same  Jacobian  may  be 
used  for  more  than  one  iteration.  In  fact,  for  weakly  nonUnear  circuits  at  low  power,  the  Jacobian 
may  need  to  be  calculated  only  once.  For  more  strongly  nonlinear  circuits,  some  criterion  may  be 
determined  for  how  often  the  Jacobian  is  updated.  TRANSIM  uses  the  same  Jacobian  for  each 
iteration  until  the  error  has  been  improved  by  less  than  5  percent  since  the  last  iterate.  At  that 
point,  the  Jacobian  is  recalculated. 

Another  method  for  reducing  run  time  is  the  blodc  Jacobian  method.  For  many  nonUnear 
circuits,  the  Jacobian  may  be  arranged  so  that  most  of  the  nonzero  elements  are  located  in  blocks 
on  the  diagonal.  Instead  of  solving  one  set  of  N  equations,  one  may  solve  M  sets  of  P  equations, 
where  M  is  the  number  of  diagonal  blocks  and  MP  =  N.  For  this  method,  (1)  becomes: 

•+ixfc  =’■  Xfc  -  aJ*n’x)  •  f(’x)  1:  =  0, 1,  •  •  • ,  M.  (3) 
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Figure  4:  Schematic  of  the  distributed  amplifier. 


Since  the  number  of  operations  required  to  solve  N  equations  is  proportionail  to  iV^/3,  it  is  clear 
that  much  computation  time  can  be  saved. 

The  problem  of  memory  usage  is  most  effectively  addressed  by  the  use  of  sparse  matrices.  For 
mzmy  problems,  the  Jacobian  is  quite  sparse,  i.e.,  there  are  many  entries  in  the  matrix  which  are 
zero.  Rizzoli  et  al.  describe  a  method  for  assuming  a  pattern  of  zero  entries  in  the  Jacobian  and 
using  specialized  matrix  solvers,  but  the  entire  matrix  still  must  be  allocated.  A  sparse  matrix 
package  such  as  the  one  developed  by  Kimdert  and  implemented  in  SPICE3  dynamically  allocates 
memory  to  relevant  entries  only.  Instead  of  storing  the  matrix  as  a  two-dimensional  array,  a  linked 
list  which  stores  entry  values  and  matrix  indices  is  used  to  store  the  same  information.  Additional 
code  is  necessary  to  manipulate  the  matrix  and  to  solve  equations.  Not  only  do  sparse  matrices 
save  memory,  but  LU  decomposition  and  matrix  equation  solving  proceed  more  quickly  due  to  the 
elimination  of  multiply  by  zero  operations. 

All  of  the  methods  discussed  above  have  been  included  in  TRANSIM  and  are  explored  in  several 
different  permutations  in  the  following  examples.  In  all  cases,  the  simulations  were  performed  on 
a  DECstation  5000. 

A. 2.1  Single  Tone  Example 

To  demonstrate  the  effectiveness  of  the  different  techniques  discussed  above,  we  consider  the 
distributed  amplifier  shown  in  Fig.  4.  First  consider  Newton’s  method  without  any  enhancements, 
i.e.,  without  sparse  matrix  routines  or  the  block  Jacobian  and  chord  methods.  The  most  time 
consuming  p2irt  of  the  iterative  process  is  building  the  Jacobian  and  finding  J"^(’x).  The  unknown 
phasor  quantities  to  be  fotmd  are  the  MESFET  terminal  voltages  and  the  current  through  each 
voltage  source  for  a  total  of  9  unknowns.  A  10  dBm  signal  at  4Ghz  is  applied  at  the  input  and 
we  consider  DC  and  7  hsirmonics  in  the  solution.  Because  the  real  emd  imaginary  parts  of  each 
quantity  are  considered  separately,  the  Jacobian  is  of  dimension  (2  •  8  ■  9)  x  (2  •  8  •  9).  The  conversion 
of  a  vector  of  time  domain  derivatives  to  the  frequency  domain  requires  a  matrix  multiplication 
followed  by  an  NFFT.  For  each  MESFET,  there  are  four  such  vectors,  yielding  twelve  matrix 
multiplications  and  transforms. 

First  the  single  tone  aincilysis  was  performed  without  any  enhancements.  Next,  the  block 
Newton  method  was  used  for  up  to  8  blocks.  The  chord  method  was  then  used  both  by  itself  and 
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Table  1:  Run  times  in  seconds  for  single  tone  analysis  of  the  distributed  amplifier. 


i 


No  Chord 

Chord 

Number  of  Blocks 

Dense 

Sparse 

Dense 

Sparse 

1 

47.6 

21.4 

24.7 

2 

35.9 

12.2 

29.1 

7.3 

4 

35.1 

9.4 

19.4 

8 

43.4 

10.9 

17.5 

Table  2:  Memory  usage  in  megabytes  for  single  tone  analysis  of  the  distributed  amplifier. 


No  Chord 

Chord 

Number  of  Blocks 

Dense 

Sparse 

Dense 

Sparse 

1 

3.9 

2.3 

3.8 

2.1 

2 

4.4 

2.3 

4.3 

2.1 

4  1 

4.1 

2.r 

3.9 

2.1 

8 

4.0 

2.3 

3.8 

2.1 

with  the  block  method.  Resulting  run  times  and  memory  usage  are  shown  in  Table  1  and  Table  2. 
Run  times  given  include  parsing  the  netlist,  initiaUzing  data  structures,  and  printing  output  as 
well  as  the  time  required  for  analysis. 

Sparse  matrix  routines  were  then  used  for  single  tone  analysis.  Both  the  block  Newton  and 
chord  methods  were  also  used  with  sparse  matrices  as  described  above.  Resulting  run  times  and 
memory  usage  are  shown  in  Table  1  and  Table  2. 

The  results  from,  each  of  the  permutations  of  solution  methods  are  the  same. 


B  Transient  analysis  of  quasioptical  systems 

This  work  is  not  complete  due  to  inadequate  modeling  of  the  active  devices  and  poor  inte^ation 
of  the  various  tools  that  have  been  developed  to  model  the  quasi-optical  system.  What  is  pre¬ 
sented  here  is  prelimina^  results.  We  expect  to  complete  the  tool  integration  m  the  subsequent 

commercialization  activities.  i.  •  i?*  i 

The  first  target  was  simulation  of  the  quasi-optical  open-cavity  oscillator  shown  in  figure  1. 

IMPATT  diodes  were  used  as  the  active  elements  and  hence  the  work  described  in  Section  2.6  on 
IMPATT  diode  modeling.  Measured  and  simulated  driving  point  impedances  are  shown  in  Figure 

5.  The  modeling  strategy  is  shown  in  Figure  6.  .  x.-  n  •  *i. 

With  just  one  diode  the  spectrum  of  the  oscillating  signal  is  shown  in  Figure  7(a)  m  the 

cavity  (soHd  line)  and  in  free  space  (dashed  line).  Also  shown  at  the  top  of  the  J! 

the  measured  driving  point  impedance  showing  a  number  of  mode  resonan^.  The  IMPATT 
diode  was  operated  in  pulsed  mode  with  the  measured  characteristics  shown  in  Figure  7(b)  were 
obtained.  The  speckled  diagram  is  the  RF  envelope  while  the  larger  pulse  is  the  ProfUe  of  the 
current  waveform.  Clearly  is  takes  considerable  time  for  osciUations  to  be  estabhshed.  The  best 
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Figure  5:  Driving  point  impedance  of  a  patch  antenna  in  a  quasi-optical  cavity;  solid  line  is 
cadculated;  dashed  line  is  measured. 


Figure  6:  Quasi-optical  system  modeling. 


simulations  to  date  are  shown  in  Figure  8  showing  that  the  correct  form  of  the  response  was 
obtadned  but  quantitative  agreement  is  lacking  largely  due  to  the  tool  integration  problems  cited 
above. 


C  Large  Signal  IMP  ATT  Modeling 

C.l  Large  Signal  IMPATT  Operation 

•  Oscillation  is  possible  with  IMPATT  diodes  because  of  a  negative  diiferential  resistcince 
provided  by  a  phase  difference  between  voltage  aind  current. 

•  For  a  reverse-biased  IMPATT  operating  under  large-signad  conditions,  there  are  two  delay 
mechanisms  which  provide  the  phase  difference. 
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Figure  7:  Measured  spectrum  and  transient  response  of  a  pulsed  IMPATT  diode  power  combiner 
in  a  Fabry  Perot  cavity:  (a)  spectrum,  solid  line  no-cavity,  dashed  line  in-cavity;  (b)  transient 


response. 
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Figure  8:  Simulated  transient  response  of  an  IMPATT  diode  power  combiner  in  a  Fabry  Perot 
cavity  using  a  simple  canonical  model  of  the  cavity. 


The  avalanche  process 

The  time  required  for  injected  charge  to  traverse  the  drift  region 

•  For  numerical  simplicity,  several  assumptions  must  be  made  for  large  signal  analysis: 

The  drift  velocity  is  the  same  for  holes  and  electrons,  and  that  velocity  is  equal  to  the 
saturation  velocity. 

Ionization  rates  are  the  same  for  holes  and  electrons. 

The  ionization  coefficient  a  can  be  written  «is 

oi{E)  = 

•  With  the  assumption  of  equal  hole  and  electron  velocities  and  ionization  rates,  the  continuity 
equations  lead  to  the  Read  equation: 

•  For  simplicity  (at  least  for  now),  assume  linear  capacitance  for  both  the  avalanche  and  drift 
regions. 

•  Only  the  ionization  integral  has  a  voltage  dependence  which  must  be  solved  in  the  time 
domaiin,  so  we  can  solve  everything  else  in  the  frequency  domain. 

Q 

•  Since  —  translates  to  multiplication  by  ju  in  the  frequency  domain,  we  can  write  the 

.  ^ 

equations  as  shown  below. 


Frequency  Domaun  Equations: 

(5) 

(6) 
(7) 


jl^kila)k  =  —  m  adx  -  (t'a  +  /o)|  + 

;wjfe(/rf)fc  =  -  {h)k) 

Td 


C.2  Goeller’s  IMPATT  device  model 

Goeller’s  model,  shown  in  Fig.  9  is  a  popular  model  for  am  IMPATT  diode  and  was  initially 
investigated  for  incorporation  in  to  the  Hainnonic  Badance  procedure.  It  was  amticipated  that  this 
would  be  a  small  drain  on  the  project.  However  this  models  turns  out  to  be  very  difficult  to 
incorporate. 

Time  Domain  Device  equations: 
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Figure  9:  Goeller’s  IMP  ATT  device  model 


^  (*load  *a) 

^Vd  =  ^{iio^d-id) 


d.  2,. 


Where 


©(E)  = 


and  the  electric  field,  E{x)  is  given  by 


R(r)  =  {  +  Va  +  Vrf  -  X<6 

1  +  Va  +'«■;  -  -  6)  X>6 


(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 


This  mathematical  model  is  difficult  to  implement  principally  because  the  time  derivative  of 
the  avalanche  current  depend  on  the  avalanche  current.  This  equation  must  be  satisfied  at  eeich 
simulation  iteration.  The  the  following  section  we  discuss  enhancements  to  the  Goeller  model  to 
mcike  it  more  suitable  for  Harmonic  Balance  analysis. 

C.3  Solving  the  Frequency  Domain  Avalanche  Current  Equation 

1  ail  -  1  j  (i.(l)  +  /o)|^  +  /,<M  (15) 
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Let  b{t)  =  /  adx  -  1  and  y{t)  =  b{t){ia{t)  +  lo)- 
Jo  ,  - 

Then  the  avalanche  current  equation  becomes; 

Mhk  =  -  (16) 

To. 

If  we  know  b{t)  and  ia{t)  +  !o  in  the  time  domain,  then  we  can  use  Chang’s  Arithmetic  Operator 
Method  (AOM)  to  find  the  spectral  content  of  y{t).  This  is  written 

y  =  Tbi  (17) 

where  Tb  is  the  “spectral  transform  matrix”  of  b{t)  and  bold  type  indicates  spectral  quantities. 

The  avalanche  current  equation  then  can  be  solved  for  AT  +  1  frequency  components  as  a  set 
oi  K  +1  linear  equations. 


0 

;wi/ai 

_ 

Boo 

Bio 

Box 

Bn 

•  •  •  Bok 
■  •  •  Bik 

foO  +  A) 
/«! 

+ 

1 

... 

ji^klaK  . 

.  Bkq 

Bki 

.  •  •  Bkk  . 

.  hK 

[o  J 

So  for  each  iteration  of  the  harmonic  balance  simulation,  we  must  solve  a  (AT  +  1)  x  (Z^"  +  1) 
system  of  equations. 

C.4  Harmonic  Balance  implementation  of  Goeller’s  IMPATT  device 
model 

•  For  simplicity  (at  least  for  now),  assume  linear  capacitance  for  both  the  avalanche  and  drift 
regions. 

•  Only  the  ionization  integral  has  a  voltage  dependence  which  must  be  solved  in  the  time 
domain,  so  we  can  solve  everything  else  in  the  frequency  domain. 

•  Since  —  translates  to  multiplication  by  ju  in  the  frequency  domain,  we  can  write  the 
equations  as  shown  below. 


Frequency  Dom2un  Equations; 


i^k{Ia)k 

j^k{Id)k 


—  l)  (t.  +  /o)|  + 

f  ({/.)*  -  (/.)») 


(19) 

(20) 
(21) 
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Figure  10:  Harmonic  balance  compatible  Goeller’s  IMPATT  device  model 


The  result  is  an  harmonic  balance  compatible  Goeller’s  IMPATT  device  model  as  shown  in 
FigurelO. 


C.5  Including  Effective  Resbtance  in  the  Frequency  Domain  Avalanche 
Current  Equation 


jUklak  =  -  1  j  (i.(i)  +  /o)|^  +  1.6(0) 


(22) 


To  solve  the  above  equation  when  the  resistance  due  to  thermal  elFects  and  space  charge  is 
included,  am  iterative  approach  must  be  used. 

At  each  harmonic  balance  iteration,  we  have  the  spectral  content  of  voltage  Va  from  the  vector 
of  imknowns  ’x. 


•  Guess  the  values  of  u,(t),  where  =  ««(<)  -  (*o(i)  +  Io)R- 

•  Using  the  v*(t)  values,  find  z„(t)  +  lo  from  the  procedure  discussed  previously. 

•  Use  these  values  in  the  equation  =  Ux(t)  +  *<>  obtain  an  error  function. 


•  Minimize  the  error  function. 


C.6  Analytical  Solution  of  the  Ionization  Integral 


We  can  rewrite  the  expression  for  electric  field: 


E{x) 


A  —  Bx  X  <b 
C  —  Dx  X  >  b 


(23) 
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where 


+Va+Vd 

B  =  ^ 

C  =  A 

D  =  3lh 

From  Maple,  the  integral 

e-/3i/.B(»)^j.  (25) 

is  evaluated  as 

j  ((-A  +  Bl)(^  +  ^iE,(hl(A  -  Bl))  +  Aeft/-'  -  A£:,(/9,M)) 

+  ^  ((i>  -  CVici^  -  0,E,{M(D  -  Cl)))  (26) 

where  Ei  is  the  exponenti<d  integral. 


C.7  Parameter  Extraction  of  the  IMPATT  Model 


To  use  the  equations  in  simulation,  we  need  VeJues  for  ai  and  j3i  which  predict  currents  and  voltages 
which  match  measured  data.  The  process  of  extracting  these  values  requires  the  formulation  of 
some  error  function  and  the  use  of  an  optimization  method  to  minimize  this  function.  A  hybrid 
optimization  method  which  uses  tree  annealing  to  find  “valleys”  of  the  error  function  and  gradient 
descent  to  explore  these  valleys  quickly  is  used  to  find  ai  ^lnd  /3i. 

The  first  attempt  to  find  these  parameters  involved  using  the  Read  equation  at  DC  to  write 
current  as  a  function  of  voltage 


where  the  measured  current  is  Iq  and  the  measured  voltage  is  present  in  E(x).  Unfortunately,  only 
one  measured  (Vo,  Iq)  pair  at  a  time  could  be  matched.  This  is  because  the  equation  as  written 
above  does  not  take  into  account  thermal  and  space  charge  resistances  which  are  present  in  the 
measurements. 

To  account  for  the  resistance  present  in  the  measurements,  the  electric  field  can  be  rewritten 


E(x)  =  Kiy/Vo  -  IqR  -  K^x 

From  the  Read  equation  we  can  write  ai  as  a  function  of 

Is/Io  +  l 


(28) 


(29) 


And  the  error  fimction  then  becoilies  the  difference  between  calculated  ori  values  for  different 
{Iq,Vq)  pairs.  This  method  of  extraction  has  found  many  different  (ai,/3i)  pairs  which  satisfy 
more  than  one  (Vo,  Iq)  pair  at  a  time. 
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